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ANALYSIS  OF  THE  PAGE  TEST  WITH  NUISANCE  PARAMETER 
ESTIMATION  FOR  VARIOUS  SIGNAL  TYPES 

INTRODUCTION 

In  [1],  the  Page  test  with  nuisance  parameter  estimation  was  introduced  as  a  method  for 
the  detection  of  a  change  from  signal- absence  (hypothesis  H)  to  signal-presence  (hypothesis 
K)  when  the  observed  data  are  statistically  characterized  by  a  model  having  unknown  noise 
parameters  with  the  same  value  under  both  H  and  K.  The  real  Gaussian  shift  in  mean 
signal  with  unknown  variance  was  considered  in  [1]  as  an  example  illustrating  the  accuracy 
of  the  Siegmund-based  approximations  to  the  average  sample  numbers  (ASNs)  for  the  Page 
test,  the  average  delay  before  detection  (D)  and  the  average  time  between  false  alarms 
(T).  In  this  report,  the  Page  test  with  nuisance  parameter  estimation  is  applied  to  the 
following  situations:  (i)  complex  Gaussian  distributed  data  with  a  shift  in  mean  signal  and 
unknown  variance,  (ii)  exponentially  distributed  data  with  an  unknown  noise  power,  and 
(iii)  multivariate  complex  Gaussian  distributed  data  with  a  deterministic  array  signal  and  an 
unknown  interference  covariance  matrix.  In  the  following  sections,  the  detector  nonlinearity, 
ASN  analysis,  and  performance  evaluation  of  each  of  the  above  situations  is  developed  and 
presented. 

A  more  appropriate  signal  model  in  sonar  consists  of  a  signal  that  starts  and  then  stops. 
The  probability  of  detecting  such  a  finite  duration  signal  using  the  Page  test  with  nuisance 
parameter  estimation  is  approximated  by  the  Brownian  motion  and  moment  matching  ap¬ 
proximations  of  [2]  and  estimated  by  the  Poisson  mixture  method  of  [3]. 

The  application  of  the  Page  test  with  nuisance  parameter  estimation  to  nonparametric 
change  detection  is  also  discussed. 
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COMPLEX  GAUSSIAN  SHIFT  IN  MEAN  SIGNAL  WITH  UNKNOWN 

VARIANCE 


Complex  Gaussian  data  are  typically  observed  at  the  output  of  discrete  Fourier  transforms 
subject  to  Gaussian  random  process  excitation  or  after  basebanding  operations.  The  complex 
shift  in  mean  signal  results  from  sinusoidal  signals  with  unknown  phase.  In  most  realistic 
situations,  the  variance  of  the  complex  Gaussian  data  is  unknown. 


STATISTICAL  DESCRIPTION 


The  observed  data  are  assumed  to  be  distributed  according  to 

CH  (0,  (T^)  i  <  q 


Xi  ~  < 


(1) 


CA/”  (/^,  cr^)  i>q 

where  cr^  is  the  unknown  variance  of  the  observed  data,  ji  is  the  unknown  complex  mean  of 
the  signal,  and  q  is  the  unknown  starting  time  of  the  signal.  The  notation  x  ~  CM  (//,  cr^) 
indicates  that  the  random  variable  x  has  a  complex  Gaussian  distribution  with  mean  /i, 
variance  and  probability  density  function  (PDF) 


f{^)  = 


1 


■KU^ 


(2) 


for  X  G  C,  where  C  is  the  complex  plane. 


DETECTOR  NONLINEARITY 

If  the  nuisance  parameter,  in  this  case  cr^,  were  known,  the  detector  nonlinearity  used  in 
the  Page  test  would  be  either  a  locally  optimal  (LO)  or  generalized  likelihood  ratio  (GLR) 
nonlinearity.  As  the  signal  is  represented  here  by  a  complex  mean,  direct  application  of 
the  LO  nonlinearity  is  not  straightforward.  However,  the  GLR  nonlinearity  may  be  easily 
derived  by  describing  the  log-likelihood  ratio  (LLR)  as 

,/  x  ,  fK{x) 


2 


TR  10,991 


\x 


X  - 


(3) 


The  GLR  statistic  is  the  maximum  of  the  LLR  over  all  possible  values  of  the  unknown 
complex  signal  amplitude  ji.  Maximizing  equation  (3)  clearly  requires  choosing  /i  =  x  which 
leads  to  the  GLR  statistic 

|2 


X 


r(x)  =  1-^. 


(4) 


Utilization  of  this  statistic  in  the  Page  test  requires  a  bias  to  force  the  mean  negative  under 
hypothesis  H  and  substitution  of  an  estimate  of  the  nuisance  parameter.  At  this  point  it  is 
useful  to  recognize  that  T  (x)  is  the  scale  of  a  noncentral  chi-squared  random  variable  with 
two  degrees  of  freedom,  thus, 

v  =  ^~xi({),  (5) 

where  S  =  ^  =  2s  is  the  noncentrality  parameter  (s  =  Thus,  the  bias  may  be  chosen 
fb“maximize“fhe  asymptotic  Page  test  performance  for  a  design  noncentrality  parameter,  S, 
as  shown  in  [4].  Unfortunately,  upon  substitution  of  the  nuisance  parameter  estimate,  this 
leads  to  very  cumbersome  analysis.  As  described  in  [4],  an  alternative  for  weak  signals  is  the 
Dyson  bias,  which  has  the  form 

T  =  n  +  ^  (6) 

for  noncentral  chi-squared  signals  with  n  degrees  of  freedom.  Thus,  upon  use  of  the  Dyson 
bias  and  substitution  of  the  nuisance  parameter  estimate,  the  Page  test  nonlinearity  has  the 
form 

g(x,a^)  = 


where  r  =  ^,a=|/i|^isa  design  signal  power,  and  s  =  ^  is  a  design  signal- to  noise  power 
ratio  (SNR). 
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The  ASN  analysis  requires  the  mean,  moment-generating-function  (MGF)  unity  root, 
and  Siegmund  correction  terms  of  the  detector  nonlinearity  (7)  under  hypotheses  H  and  K 
conditioned  on  the  observed  nuisance  parameter  estimate  The  detector  nonlinearity  is 
a  scaled  and  shifted  noncentral  chi-squared  random  variable  under  H  and  K.  A  numerical 
routine  for  obtaining  the  MGF  unity  root  is  presented  in  Appendix  A.  The  conditional  mean 
is 


E 


\  (  -2^ 

-  2 

g  (2),  (T  ^ 

E[r] 

r 

2  +  5 


-2-  - 
r 


1  +  s 


-2, 


(9) 


where  here  and  in  the  following,  hypothesis  H  is  obtained  by  setting  5  =  5  =  0.  The 
Siegmund  correction  terms  conditioned  on  are 


P±  =  —  , 
r 


(10) 


where  are,  as  described  in  [1],  the  Siegmund  correction  terms  for  a  noncentral  chi-squared 
random  variable  with  two  degrees  of  freedom  and  noncentrality  parameter  6.  As  found  in  [1], 
the  correction  terms  when  no  signal  is  present  are  p°  =  2  and  p°  =  —  |. 


EXPECTATIONS  OVER  NUISANCE  PARAMETER  ESTIMATE 

As  seen  in  [1],  the  Siegmund-based  approximation  to  the  ASN  of  the  Page  test  with 
nuisance  parameter  estimation  requires  the  expectation  of  three  functions  of  the  nuisance 
parameter  estimate.  They  are  described  in  this  section  for  the  complex  Gaussian  shift  in 
mean  signal.  The  nuisance  parameter  estimate  is  the  maximum  likelihood  estimate  (MLE) 
of  the  variance  over  the  auxiliary  data.  If  the  auxiliary  data  (yi, . . .  ,yj)  are  independent 
and  identically  distributed  according  to 

%  -  CA/"  (O,  cr^)  ,  (11) 

the  nuisance  parameter  estimate  has  the  form 

=  (12) 

^  J=1 
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Since 


2M 


^2  is  centrally  chi-squared  distributed  with  two  degrees  of  freedom,  2J^ 


Thus,  r  follows  a  Gamma  distribution 


r  =  ^  ~  Gamma  (j,-z]  , 


V  ’  J 


(13) 


which  has  PDF 


fir)  = 


‘■Q 

~T{J) 


1^-Jt 


(14) 


for  r  >  0. 

The  first  expectation  is 


Elis)  =.E  : 


-  — E 


2  [r-(l  +  s-|)J 


(15) 

(16) 


where,  if  r  follows  the  gamma  distribution  with  degrees  of  freedom  parameter  a  and  scale 
parameter  /?,  the  expectation  of  equation  (15)  may  be  described  by  the  general  function 


7  (7,  a,/?)  =  E 


r  —  7 


-  /,.o( 


r  —  7/  r (a) 

_  r  MfL,. 

~  7r.or'  '' 


dr 


(17) 

(18) 


jE  (a)  (r  —  7) 

with  the  notation  f  representing  the  principal  value  integral.  By  moving  the  principal  value 
integral  of  equation  (18)  from  the  real  axis  upward  into  the  complex  plane,  7 (7, a,/?)  may 
be  described  as 


l(~l,a,0)  =  f  / (z) dz  +  wRes {/ (z) ; 2  =  7} 


(19) 


where 


fiz)  = 


/  \  Ct  z 

(f)  °  ^ 

r(a)(^-7)  ’ 


(20) 
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and  C  is  any  contour  in  the  first  quadrant  from  the  origin  to  z  =  oo.  The  residue  of  /  (z)  at 
z  =  7  is 


Res  {/  (z) ;  z  =  7}  = 


The  contour  C  may  be  described  as  the  limit  as  j4  00  of  U  Ca  where  Cq^a  is  the 
integral  from  the  origin  to  the  point  and  Ca  is  the  integral  along  the  arc  of  radius  A 
from  angle  d  to  angle  zero.  Due  to  the  exponential  term  in  /  (z),  this  latter  integral  goes  to 
zero  as  ^  ^  00.  Thus,  letting  z  =  (dz  =  eP^da)  results  in 


s=7 


-V 

0) 


r(a) 


(21) 


[  f{z)dz  =  f  f{z)dz 

JC  JCb,oo 

=  /  (ae^'^)  da  ,  (22) 

_  ^ 

which,  when  9  is  appropriately  chosen,  may  be  easily  performed  numerically.  Due  to  the  e  ^ 
term  in  /  (z),  the  substitution  z  =  results  in  the  sinusoidal  term  e  ^  ^  .  This  function 

causes  problems  in  numerical  integration  routines  when  jd  is  small.  Thus,  it  is  suggested  that 
9  be  chosen  proportional  to  fd  to  keep  the  frequency  of  the  sinusoidal  terms  approximately 
constant  as  fd  varies.  Choosing  9  =  was  observed  to  work  well  for  (d  between  0.1  (J  =  10) 
and  0.01  (J  =  100). 

It  should  be  noted  that  I  (7,  a,  (d)  may  be  computed  recursively  if  a.  is  an  integer  with  an 
initial  starting  value  requiring  evaluation  of  the  exponential  integral.  However,  the  recursion 
is  unstable  for  large  values  of  a  (i.e.,  a  =  J  >  30). 

The  second  expectation  is 


Et{s)  =  E 


P± 


[E[g\a^ 


Pi 


E 


Now,  observe  that  equation  (17)  may  be  written  as 

•00  ^ 

r  (a)  (r  —  7) 


/oo 

- 

r=0 


dr 


(23) 
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roo  0  foo 

X=or(Q)/?“'^''^  V  r=0 


rp(^  ^0  0 


=  1  +  7E 


r  (q)  /?“ 

1 


r  (a)  /3“  (r  —  7) 


dr 


r  —  7 


(24) 


where  the  first  integral  in  the  second  line  is  recognized  as  the  integral  of  the  PDF  of  r.  Thus, 

. 


E 


1 

r  —  7 


7 


(25) 


and 


£^2"^  (5)  =  - 


Pi 


2(1 


+  5 


(26) 


The  third  expectation  has  the  form 


E^{s)  =  E 

=  piE 


tsP- 


[1  _  g4,(/i+p+-p-)  J 


[1  _gt,(/i+r-l(p^-pl))J 

5  f°°  r~^ts 

~  Jr=0  I  _  (*+»•- ^(p+-p1)) 


/  {r)  dr  , 


(27) 


where  ts  is  the  MGF  unity  root  of  the  Page  test  update  when  5  =  2s  is  the  noncentrality 
parameter  and  the  scale  and  bias  terms  are  as  found  in  equation  (8).  This  expectation  is 
easily  performed  by  standard  numerical  integration  with  tg  obtained  as  a  function  of  the 
variables  (s,  s,  r)  using  the  method  described  in  Appendix  A. 


AVERAGE  SAMPLE  NUMBER  AND 
DETECTION  PROBABILITY  ANALYSIS 

Combining  equations  (16),  (26),  and  (27)  as  described  in  [1]  results  in  the  Siegmund-based 
approximation 


(5) 


hEi  (s)  +  E2  (s)  —  E2  (s)  + 


E2  is) 

Ez{s) 


~h 

~Y 


I 


2iE^{s) 


(28) 


where  7  =  1  +  s  —  The  average  time  between  false  alarms  is  obtained  when  5  =  s  =  0. 
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When  the  nuisance  parameter  is  perfectly  estimated  (i.e.,  J  — *•  oo),  r  — >•  1  and  the  Page 
test  with  nuisance  parameter  estimation  results  in  the  standard  Page  test.  In  this  case,  the 
Siegmund-based  approximation  to  the  ASN  is 


1  +  tshg  — 
2ts  (^s  — 


(29) 


where  hg  =  h  +  -  pt  and  tg  is,  as  before,  the  MGF  unity  root  where  r  =  1  in  the  scale 

and  bias  terms  of  equation  (8). 

The  above  Siegmund-based  approximations  to  the  average  time  between  false  alarms 
(5  =  s  =  0)  and  the  average  delay  before  detection  (6  =  2s  >  0)  are  used  to  evaluate  the 
performance  of  the  Page  test  with  nuisance  parameter  estimation  for  the  complex  Gaussian 
shift  in  mean  signal.  Figure  1  contains  the  Page  test  operating  characteristic  (OC)  curve  for 
s  =  0  dB,  a  design  SNR  s  =  0  dB,  and  for  auxiliary  data  amounts  J  =  10,  20,  30,  50,  100, 
and  oo  (perfect  estimation  of  the  nuisance  parameter).  Curves  of  the  average  time  between 
false  alarms  as  a  function  of  the  threshold  for  the  above  situation  are  found  in  figure  2.  The 
average  delay  before  detection  as  a  function  of  SNR  is  found  in  figure  3  where  the  design 
SNR  and  average  time  between  false  alarms  are  held  constant  at  s  =  0  dB  and  T  —  10^, 
respectively.  The  signal  level  that  produces  a  zero  mean  in  the  Page  test  update,  and 
therefore  a  very  large  average  delay  before  detection,  is  considered  the  minimum  detectable 
signal  level  (MDL).  In  this  situation  the  MDL  increases  as  J  decreases  from  approximately 
—3  dB  for  J  —  oo. 


Simulations  are  performed  to  validate  the  approximations  to  T  and  D  and  to  evaluate 
the  performance  of  existing  methods  for  approximating  the  probability  of  detecting  a  finite- 
duration  signal  using  the  Page  test  with  nuisance  parameter  estimation.  Figures  4  and  5 
contain  the  Siegmund-  and  Wald-based  approximations  to  T  and  D  along  with  that  estimated 
from  simulation.  As  seen  in  [1],  the  Wald-based  approximations  may  be  obtained  from  Ei  (s) 
and  from  (s)  with  a  slight  modification.  As  the  Siegmund-based  approximations  are  more 
accurate  than  the  Wald-based  ones,  this  modification  is  left  to  those  interested  in  analyzing 
the  Wald-based  approximations.  In  all  cases,  the  design  SNR  is  0  dB  and  J  =  30  samples 
of  auxiliary  data  are  used  to  estimate  the  nuisance  parameter  with  a  60  sample  buffer.  The 
thresholds  for  the  results  shown  in  figure  5  are  set  according  to  T  =  10^.  The  discrepancy 
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between  the  Siegmund-based  approximation  to  D  and  the  simulations  at  lower  SNR  is  the 
result  of  an  inadequate  buffer  size.  In  this  situation,  the  nuisance  parameter  estimate  is 
corrupted  by  data  containing  signal,  particularly  for  SNRs  below  1.5  dB. 

The  most  common  use  for  the  Page  test  is  for  the  detection  of  a  signal  onset  where,  once 
the  signal  starts,  it  does  not  stop.  The  Page  test  has  also  been  suggested  for  the  detection 
of  a  temporary  change;  that  is,  for  detection  of  finite-duration  signals  [2],  [5],  [6].  As  de¬ 
scribed  in  [2]  and  [3],  the  probability  of  detecting  such  signals  may  be  approximated  through 
simulation  based  or  analytical  methods.  Figure  6  contains  the  probability  of  detection  ap¬ 
proximated  by  the  sample  cumulative  distribution  function  (CDF)  using  2000  samples  of 
the  stopping  time  random  variable,  the  Brownian  motionT  and  moment  matching  methods 
of  [2]  and  the  Poisson  mixture  method  of  [3]  using  500  samples  for  T  =  10^,  a  design  SNR 
of  0  dB,  and  0  and  5  dB  signal  strengths.  Similar  to  the  results  of  [3],  the  Poisson  mixture 
method  provided  the  best  performance  as  is  also  indicated  by  the  total  squared  error  between 
the  latter  three  methods  and  the  sample  CDF  found  in  figure  7  as  a  function  of  SNR.  The 
decay  in  the  Brownian  motion  and  moment  matching  approximations  may  be  attributed  to 
the  corruption  of  the  auxiliary  data  by  signal  presence  for  lower  SNR  and  to  the  Page  test 
update  being  non-Gaussian  and  a  dependent  process  over  time. 

The  analytical  methods  of  [2]  require  the  unconditional  mean  and  variance  of  the  Page 
test  update.  It  can  be  shown  that  for  the  complex  Gaussian  signal  with  nuisance  parameter 
estimation, 

EM  =  +  (3°) 


and 


Var  [g]  = 
where,  if  r  Gamma  j)) 


4J2 


(J-l)(J-2) 


1  +  25  + 


(l  +  5  -  f) 
J- 1 


2  ^ 


J 

J-1 


and 


(31) 


(33) 
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have  been  used. 


logto(Average  Time  Between  False  Alarms) 


Figure  1.  Average  Delay  Before  Detection  Vs.  Average  Time  Between  False 
Alarms  for  0  dB  SNR  and  Varying  Amounts  of  Auxiliary  Data 
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Threshold 


Figure  2.  Average  Time  Between  False  Alarms  Vs.  Threshold  for  Varying 

Amounts  of  Auxiliary  Data 


SNR  (dB) 


Figure  3.  Average  Delay  Before  Detection  Vs.  SNR  for  T  =  10^  Samples  and 

Varying  Amounts  of  Auxiliary  Data 
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Figure  4.  Simulation  and  Siegmund-  and  Wald-Based  Approximations  to 
Average  Time  Between  False  Alarms  Vs.  Threshold  for  J  =  30 


Figure  5.  Simulation  and  Siegmund-  and  Wald-Based  Approximations  to 
Average  Delay  Before  Detection  Vs.  SNR  for  J  =  30  and  T  =  10^  Samples 
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Figure  6.  Approximations  to  the  Probability  of  Detection  Vs.  Signal  Duration 
for  0  dB  and  5  dB  SNR,  J  =  30,  and  f  =  10^  Samples 
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Figure  7.  Total  Error  Between  Sample  CDF  and  Probability  of  Detection 
Approximations  Vs.  SNR  for  J  =  30  and  T  =  10^  Samples 
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EXPONENTIAL  SIGNAL  WITH  UNKNOWN  NOISE  POWER 


Exponential  data  are  typically  observed  as  the  magnitude  squared  output  of  discrete  Fourier 
transforms  subject  to  zero-mean  Gaussian  random  process  excitation.  When  a  Gaussian 
signal  is  present  in  the  time  series  data,  the  mean  of  the  exponential  data  increases  by  some 
unknown  amount  6  from  the  nominal  noise  power  level  A.  In  most  realistic  situations,  the 
noise  power  A  of  the  exponential  data  is  unknown. 


STATISTICAL  DESCRIPTION 


The  observed  data  are  assumed  to  be  distributed  according  to 


Xi  ~  < 


Exp (A) 
Exp  (A  +  6) 


i  <  q 
i  >  q 


(34) 


where  q  is  the  unknown  starting  time  of  the  signal.  The  notation  x  ~  Exp  (A)  indicates  that 
the  random  variable  x  has  an  exponential  distribution  with  mean  A  and  PDF 


f  (x)  = 


(35) 


for  X  >  0. 


DETECTOR  NONLINEARITY 

If  the  nuisance  parameter,  in  this  case  A,  were  known,  the  detector  nonlinearity  used 
in  the  Page  test  would  be  either  the  LLR  or  the  LO  nonlinearity.  Here,  the  LLR  and  LO 
nonlinearities  result  in  the  same  form. 


(x)  =  log 


fKjx) 

fnix) 

9 


A  (^  -|-  A) 


X 


log  1  +  - 


(36) 


14 


TR  10,991 


Assuming  that  6  is  unknown  results  in  scaling  equation  (1567T)y  ^  andchoosmg  a  design 
signal  strength  $  for  the  bias  term, 


l{x) 


X 

X 


(37) 


This  is  equivalent  to  using  the  asypmtotically  optimal  bias  [4]  on  the  nonlinearity  f .  As 
in  the  previous  section,  the  Dyson  bias  results  in  simpler  analysis.  The  Dyson  bias  for  the 
above  nonlinearity  has  the  form 


r 


(38) 


where  s  =  |  is  a  design  SNR  and 


r 


A 

A 


(39) 


is  the  normalized  estimate  of  the  background  noise  power. 

Thus,  the  Page  test  nonlinearity  with  the  Dyson  bias  and  substitution  of  the  nuisance 
parameter  estimate  has  the  form 


.Ra)  =  ?-l-± 

r  2r 


(40) 

(41) 


where  y  =  |  ~  Exp  (1  +  s)  where  5  =  |  is  the  actual  SNR. 

The  ASN  analysis  requires  the  mean,  variance,  MGF  unity  root,  and  Siegmund  correc¬ 
tion  terms  of  the  detector  nonlinearity  (40)  under  hypotheses  H  and  K  conditioned  on  the 
observed  nuisance  parameter  estimate  A.  The  detector  nonlinearity  is  a  scaled  and  shifted 
exponential  random  variable  under  H  and  K.  A  numerical  routine  for  obtaining  the  MGF 
unity  root  is  presented  in  Appendix  B.  The  conditional  mean  is 

E[R] 


E  [c?  |A  = 


r 

1  +  5 


^  2r 
T  2r 

^  5 


(42) 
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where  here  and  in  the  following,  hypothesis  H  is  obtained  by  setting  s  =  0.  The  conditional 
variance  is 


Var 


(z,  a) 


ivarll/] 

(1  +  T 

r2 


The  Siegmund  correction  terms  conditioned  on  A  are 

p±  =  (1  +  s)  —  , 
r 

where  =  1  and  p°  =  —  |  are  as  described  in  [1]. 


(43) 


(44) 


EXPECTATIONS  OVER  NUISANCE  PARAMETER  ESTIMATE 


The  three  expectations  over  the  nuisance  parameter  estimate  required  for  the  Siegmund- 
based  ASN  approximation  are  described  in  this  section  for  the  exponential  signal.  The 
nuisance  parameter  estimate  is  the  MLE  of  the  noise  power  over  the  auxiliary  data.  If  the 
auxiliary  data  (yi, . . .  ,yj)  are  independent  and  identically  distributed  according  to 


yj  ~  Exp  (A) 


(45) 


the  nuisance  parameter  estimate  has  the  form, 

1  ^ 

which,  when  scaled  by  ^  is  easily  shown  to  follow  a  gamma  distribution 

r  =  ^  ~  Gamma  . 

Using  the  result  of  equation  (18),  the  first  expectation  is 

E,{s)  =  e[e-‘[s|a]' 

=  -E 


[r-  (l  +  s-|)J 


(46) 


(47) 


(48) 
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Using  the  result  of  equation  (25),  the  second  expectation  is 

P± 


E^{s)  =  E 


[e  [g\X 
(1  +  S)P±E 


1 


_r-  (l  +  s-f) 


=  -  (1  +  s) 

The  third  expectation  has  the  form 
E,{s)  =  E 


/(l-hs-  f,  J,  7 

1 

(l  +  s-  f) 

tsp- 


_  Qts{h+p+-p-) 

1  +  S 


1  +  S 


-E 


r  H. 


_  Q^s^h+r  ^(H-s)(p^-p2,))  J 

X=0  I  _  g*s(^+^(l+*)) 


-/  (r)  dr 


(49) 


(50) 


where  tg  is  the  MGF  unity  root  of  the  Page  test  update  when  the  scale  and  bias  terms  are 
as  found  in  equation  (41).  This  expectation  is  easily  performed  by  standard  numerical  inte¬ 
gration  with  ts  obtained  as  a  function  of  the  parameters  (s,  s,  r)  using  the  method  described 
in  Appendix  B. 


AVERAGE  SAMPLE  NUMBER  AND 
DETECTION  PROBABILITY  ANALYSIS 


Combining  equations  (48),  (49),  and  (50)  as  described  in  [1]  results  in  the  Siegmund-based 
approximation 


As  (s) 


hE\  (s)  -l-  E2  (s)  —  E2  (s)  + 


e;  (s) 

E,(s) 


(51) 


where  7  =  1  -H  s  —  The  average  time  between  false  alarms  is  obtained  when  s  =  0. 

When  the  nuisance  parameter  is  perfectly  estimated  (i.e.,  J  00),  r  —>■  1  and  the  Page 
test  with  nuisance  parameter  estimation  results  in  the  standard  Page  test.  In  this  case,  the 
Siegmund-based  approximation  to  the  ASN  results  is 


1  +  tshs  - 

ts  {s  -  I) 


(52) 
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where 

hg  =  h  +  p_|_  —  p_ 

—  ^  +  -  (1  +  s)  ,  (53) 

and  tg  is,  as  before,  the  MGF  unity  root  where  r  =  1  in  the  scale  and  bias  terms  of  equa¬ 
tion  (41). 

The  above  Siegmund-based  approximations  to  the  average  time  between  false  alarms 
(s  =  0)  and  the  average  delay  before  detection  (s  >  0)  are  used  to  evaluate  the  performance 
of  the  Page  test  with  nuisance  parameter  estimation  for  the  exponential  signal.  Figure  8 
contains  the  Page  test  OC  curve  for  s  =  0  dB,  a  design  SNR  s  =  0  dB,  and  for  auxiliary 
data  amounts  J  =  10,  20,  30,  50,  100,  and  oo  (perfect  estimation  of  the  nuisance  parameter). 
Curves  of  the  average  time  between  false  alarms  as  a  function  of  the  threshold  for  the  above 
situation  are  found  in  figure  9.  The  average  delay  before  detection  as  a  function  of  SNR  is 
found  in  figure  10  where  the  design  SNR  and  average  time  between  false  alarms  are  held 
constant  at  s  =  0  dB  and  T  =  10^  respectively.  Note  that  in  this  situation  there  is  a  MDL 
that  increases  as  J  decreases  (note  that  the  MDL  for  J  =  oo  is  approximately  —3  dB). 

Simulations  are  performed  to  validate  the  approximations  to  f  and  D  and  to  evaluate 
the  performance  of  existing  methods  for  approximating  the  probability  of  detecting  a  finite- 
duration  signal  using  the  Page  test  with  nuisance  parameter  estimation.  Figures  11  and  12 
contain  the  Siegmund-  and  Wald-based  approximations  to  T  and  D  along  with  that  estimated 
from  simulation.  In  all  cases,  the  design  SNR  is  0  dB  and  J  =  30  samples  of  auxiliary  data  are 
used  to  estimate  the  nuisance  parameter  with  a  buffer  of  equivalent  length.  The  thresholds 
for  the  results  shown  in  figure  12  are  set  according  to  T  =  10^.  The  discrepancy  between 
the  Siegmund-based  approximations  and  the  simulations  at  lower  SNR  is  the  result  of  an 
inadequate  buffer  size.  In  this  situation,  the  nuisance  parameter  estimate  is  often  corrupted 
by  data  containing  signal,  particularly  for  SNRs  below  1  dB.  Unfortunately,  the  solution 
to  this  problem  suggested  in  [1],  requiring  estimation  of  the  nuisance  parameter  under  the 
signal-present  hypothesis,  is  not  suitable  for  the  exponential  signal.  Thus,  the  buffer  zone 
separating  the  auxiliary  data  from  the  current  test  sample  must  be  carefully  designed. 

Figure  13  contains  the  probability  of  detection  approximated  by  the  sample  CDF  using 
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2000  samples  of  the  stopping  time  random  variable,  the  Brownian  motion  and  moment 
matching  methods  of  [2]  and  the  Poisson  mixture  method  of  [3]  using  500  samples  for  T  =  10^, 
a  design  SNR  of  0  dB,  and  0  and  5  dB  signal  strengths.  Similar  to  the  results  of  [3],  the 
Poisson  mixture  method  provided  the  best  performance  as  is  also  indicated  by  the  total 
squared  error  between  the  latter  three  methods  and  the  sample  CDF  found  in  figure  14  as  a 
function  of  SNR. 

The  analytical  methods  of  [2]  require  the  unconditional  mean  and  variance  of  the  Page  test 
update.  It  can  be  shown  that  for  the  exponential  signal  with  nuisance  parameter  estimation, 

EM  =  (l+.-£)(^)-l  (54) 


and 


Var  [p]  = 


(J-l)(J-2) 


,  (l  +  s-|)' 
{l  +  sf+^ - ^ 


J-  1 


where,  equations  (32)  and  (33)  have  been  used. 


(55) 
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Figure  9.  Average  Time  Between  False  Alarms  Vs.  Threshold  for  Varying 

Amounts  of  Auxiliary  Data 


Figure  10.  Average  Delay  Before  Detection  Vs.  SNR  for  T  =  10^  Samples  and 

Varying  Amounts  of  Auxiliary  Data 
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Figure  11.  Simulation  and  Siegmund-  and  Wald-Based  Approximations  to 
Average  Time  Between  False  Alarms  Vs.  Threshold  for  J  =  30 


Figure  12.  Simulation  and  Siegmund-  and  Wald-Based  Approximations  to 
Average  Delay  Before  Detection  Vs.  SNR  for  J  —  30  and  T  —  10^  Samples 
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Signal  Duration 


Figure  13.  Approximations  to  the  Probability  of  Detection  Vs.  Signal  Duration 
for  0  dB  and  5  dB  SNR,  J  =  30,  and  f  =  10^  Samples 
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Figure  14.  Total  Error  Between  Sample  CDF  and  Probability  of  Detection 
Approximations  Vs.  SNR  for  J  =  30  and  T  =  10^  Samples 
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MULTIVARIATE  COMPLEX  GAUSSIAN  DETERMINISTIC  SIGNAL 
WITH  UNKNOWN  INTERFERENCE  COVARIANCE  MATRIX 


The  multivariate  complex  Gaussian  deterministic  signal  arises  as  an  extension  of  the  complex 
Gaussian  shift  in  mean  signal  to  multiple  sensor  observations.  The  unknown  variance  extends 
to  an  unknown  covariance  matrix.  The  multivariate  deterministic  signal  model  represents 
a  sinusoidal  signal  with  unknown  phase  received  by  an  array  of  sensors  and  results  in  an 
unknown  complex  amplitude  scaling  a  known  direction  or  steering  vector  as  the  mean  of  the 
vector  observation. 


STATISTICAL  DESCRIPTION 


The  observed  data  are  assumed  to  be  distributed  according  to 


Xi  ~  < 


CUn{0,'E) 

CUn  (^d,  S) 


i  <  q 
i  >  q 


(56) 


where  S  is  the  unknown  interference  covariance  matrix  of  the  observed  data,  9  is  the  unknown 
complex  amplitude  of  the  signal,  d  is  a  known  array  steering  vector,  and  q  is  the  unknown 
starting  time  of  the  signal.  The  notation  x  ~  CMn  (^,  S)  indicates  that  the  length  N 
random  vector  x  has  a  multivariate  complex  Gaussian  distribution  with  mean  fi,  covariance 
matrix  S,  and  PDF 


/(X)  = 


TT 


(57) 


for  X  6  C^,  where  the  superscript  ^  indicates  the  complex  conjugate  and  transpose  operation. 
It  is  assumed  that  the  steering  vector  d  has  unit  length,  d-^d  =  1. 


DETECTOR  NONLINEARITY 

If  the  nuisance  parameter  S  were  known,  the  detector  nonlinearity  used  in  the  Page  test 
would  be  either  a  LO  or  GLR  nonlinearity.  As  with  the  univariate  case,  direct  application 
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of  the  LO  nonlinearity  is  not  straightforward.  However,  the  GLR  nonlinearity  may  be  easily 
derived  by  describing  the  LLR  as 

/k(x) 


=  log  >  ,  s 
fn  (x) 

=  x^S-^x  -  (x  -  (x  -  ^d) 


(x) 


(58) 

The  GLR  statistic  is  the  maximum  of  the  LLR  over  all  possible  values  of  the  unknown 
complex  signal  amplitude  6.  Maximizing  equation  (58)  clearly  requires  choosing  0  = 
which  leads  to  the  GLR  statistic 


d^S-^x 

2 

_ 

d^S-^x 

di/s-id 

^  d^S-id 

V  =  2 


d^s-i 


d^E-id  ■ 

Using  Theorem  1  of  [7]  it  is  straightforward  to  show  that  if  x  ~  CM n  (z^,  S)  then 


(59) 


w^x 


w-^Sw 

where  the  noncentrality  parameter  is 

6  =  2 


xi  («)  . 


(60) 


W^fJ, 


'w^Sw ■ 

Applying  this  to  equation  (59)  with  w  =  it  is  seen  that 

where 

6  =  2\0fd^T,-^d 

=  2s. 


(61) 


(62) 


(63) 


Thus,  the  bias  required  for  application  of  this  test  statistic  to  the  Page  test  may  be  chosen 
in  the  same  manner  as  the  univariate  case  (equation  (6)). 

Upon  substitution  of  the  nuisance  parameter  estimate,  S,  the  Page  test  nonlinearity  has 
the  form 


9 


2 


d^S-ix|^ 

d"E-id 


-  2  -  ad^t-^d 


(64) 
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W-^X, 
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'  d^E-id 
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d^E-M 
d^E-id  ’ 


(65) 


where  a  =  |(9|  is  a  design  signal  power,  s  —  ad'^E"'M  is  a  design  signal-to-interference 
power  ratio  (SIR),  and  w  =  S“M.  By  applying  equations  (60)  and  (61),  it  is  seen  that  V  is 
noncentral  chi-squared  distributed, 

|2 


V  =  2 


w^x 


w^Ew 


^2  {S)  , 


with  noncentrality  parameter 


6^2 


(^d) 
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Defining 


and 


d^E-id 

d^S-id 


(d^S-M 


(d^E-iES-id)  (d^E-id)  ’ 
the  Page  test  nonlinearity  may  be  written  as 

g(x,t)  =  --2-- 
\  /  pr  r 


(66) 


(67) 


(68) 


(69) 


(70) 


and  the  noncentrality  parameter  of  P  as  <5  =  2sp. 

The  ASN  analysis  requires  the  mean,  MGF  unity  root,  and  Siegmund  correction  terms  of 
the  detector  nonlinearity  (70)  under  hypotheses  H  and  K  conditioned  on  the  observed  nui¬ 
sance  parameter  estimate  S.  As  the  detector  nonlinearity  is  a  scaled  and  shifted  noncentral 
chi-squared  random  variable  under  H  and  K,  the  numerical  routine  presented  in  Appendix  A 
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may  be  used  for  obtaining  the  MGF  unity  root.  The  conditional  mean  is 


E  [9  (l,  S 


EM  „  5 

pr  r 

2  T  2sp 


pr 


-2-  - 
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2 

pr  L 


-2, 


where  here  and  in  the  following,  hypothesis  H  is  obtained  by  setting  6  =  s  = 
Siegmund  correction  terms  conditioned  on  S  are 


(71) 
0.  The 


P± 


4 

pr 


(72) 


where  p^.  are  as  previously  described,  except  here  8  =  2sp. 


EXPECTATIONS  OVER  NUISANCE  PARAMETER  ESTIMATE 


The  three  expectations  over  the  nuisance  parameter  estimate  required  for  the  Siegmund- 
based  ASN  approximation  are  described  in  this  section  for  the  multivariate  complex  Gaussian 
data  with  a  deterministic  array  signal.  The  nuisance  parameter  estimate  is  the  MLE  of  the 
interference  covariance  matrix  over  the  auxiliary  data.  If  the  auxiliary  data  (yi,  •  • . ,  yj)  are 
independent  and  identically  distributed  according  to 

~  CM N  (0,  S)  ,  (73) 

the  nuisance  parameter  estimate  has  the  form 

J=l 

which  is  distributed  according  to  the  scale  of  a  complex  Wishart  distributed  matrix, 

JS  ~CWiv(J,S)  ,  (75) 


where  the  notation  A  ~  CWjv  indicates  that  the  iV-by-V  random  matrix  A  has  a 

complex  Wishart  distribution  with  J  degrees  of  freedom  and  S  for  a  scale  matrix.  Utilizing 
this  and  the  results  stated  in  Appendix  A  of  [8]  it  is  seen  that  r  of  equation  (68)  is  gamma 
distributed 


Gamma 


(j-A  +  l, 


(76) 
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with  PDF 


/  W  = 


^  '  r(j-iV  +  i) 

for  r  >  0,  and  that  p  of  equation  (69)  is  beta  distributed 


p  ~  Beta  {J  —  N  +  2,  N  —  1) 


with  PDF 


for  0  <  p  <  1. 


fip)  = 


r{J  +  l)p^-^+^{l-pf- 
T{J  -  N  +  2)r{N  -1) 


Using  the  result  of  equation  (18),  the  first  expectation  is 


E,{s)  =  Ej,  E-MpiS 


where  subscripts  on  expectations  indicate  the  random  variable  over  which  the  expectation 
is  performed. 

Using  the  result  of  equation  (25),  the  second  expectation  is 
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The  third  expectation  has  the  form 
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where  tg  is  the  MGF  unity  root  of  the  Page  test  update  when  S  =  2sp  is  the  noncentrality 
parameter  and  the  scale  and  bias  terms  are  as  found  in  equation  (70).  This  expectation  is 
easily  performed  by  standard  numerical  integration  with  tg  obtained  as  a  function  of  the 

parameters  (s.s.r,  p)  using  the  method  described  in  Appendix  A. _ 

Note  that  when  A’’  =  1,  p  becomes  a  degenerate  random  variable  with  value  one  (p  — >•  1) 
and  all  of  the  results  of  this  section  simplify  to  those  of  the  univariate  case  previously 
described. 


AVERAGE  SAMPLE  NUMBER  ANALYSIS 


Combining  equations  (80),  (81),  and  (82)  as  described  in  [1]  results  in  the  Siegmund-based 
approximation 


A5  (s)  =  hEi  (s)  +  E2  (s)  —  E2  (s)  + 


Ei  is) 


2 


-E, 


E,(.s) 

/(iV-W+l.i)] 

BsWJI  2p7 


0^ 


(83) 


where  7  =  ^  +  5  —  f .  The  average  time  between  false  alarms  is  obtained  when  6  =  s  =  0. 

When  the  nuisance  parameter  is  perfectly  estimated  (i.e.,  J  — >  00),  r  1  and  p  ^  1  and 
the  Page  test  with  nuisance  parameter  estimation  results  in  the  standard  Page  test.  In  this 
case,  the  Siegmund-based  approximation  to  the  ASN  is  identical  to  that  of  equation  (29) 
with  s  = 

The  above  Siegmund-based  approximations  to  the  average  time  between  false  alarms 
(s  =  0)  and  the  average  delay  before  detection  (s  >  0)  are  used  to  evaluate  the  performance 
of  the  Page  test  with  nuisance  parameter  estimation  for  the  complex  Gaussian  shift  in  mean 
signal.  In  all  cases,  the  adaptive  dimension  is  A  =  5  and  the  design  SNR  is  0  dB.  First,  the 
approximations  are  compared  to  simulation  in  figures  15  and  16  for  J  =  30  as  a  function 
of  threshold.  The  Siegmund-based  approximations  for  varying  amounts  of  auxiliary  data 
(J  =  25,  30,  40,  50,  60,  75,  100,  and  00)  as  a  function  of  threshold  are  found  in  figures  17 
and  18. 
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Figure  18  contains  the  Page  test  OC  curve  for  s  =  0  dB  and  a  design  SNR  s  =  0  dB 
for  auxiliary  data  amounts  (perfect  estimation  of  the  nuisance  parameter).  Curves  of  the 
average  time  between  false  alarms  as  a  function  of  the  threshold  for  above  situation  are 
found  in  figure  17.  Curves  of  the  average  delay  before  detection  as  a  function  of  threshold 
are  found  in  figure  18. 
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Figure  15.  Simulation  and  Siegmund-Based  Approximation  to  Average  Time 
Between  False  Alarms  Vs.  Threshold  for  A  =  5  and  J  =  30 


Figure  16.  Simulation  and  Siegmund-Based  Approximation  to  Average  Delay 
Before  Detection  Vs.  Threshold  for  V  =  5,  J  =  30,  and  0  dB  SNR 
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Figure  17.  Average  Time  Between  False  Alarms  Vs.  Threshold  for  0  dB  SNR. 
iV  =  5,  and  Varying  Amounts  of  Auxiliary  Data 


Figure  18.  Average  Delay  Before  Detection  Vs.  Threshold  for  0  dB  SNR, 
N  —  5,  and  Varying  Amounts  of  Auxiliary  Data 
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NONPARAMETRIC  CHANGE  DETECTION 


The  nuisance  parameter  estimation  concept  of  [1]  may  be  extended  to  the  nonparametric 
detection  problem.  In  this  situation,  the  PDF  of  the  observed  data  is  unknown  under  both 
the  signal-absent  and  signal-present  hypotheses.  The  signal-present  hypothesis  is  solely  dis¬ 
tinguished  from  the  signal-absent  hypothesis  by  general  functional  relationships  between  the 
signal-absent  CDF,  F  (x),  and  the  signal- present  CDF,  G  (x).  Common  relationships  include 
a  change  in  location  parameter,  G  (x)  =  F  {x  -  6),  change  in  scale  parameter,  G{x)  =  F 
or  relationships  such  as  Lehmann  alternatives  G  (x)  =  [F  (x)T  where  ^  >  0  and  not  equal  to 
one.  Thus,  the  nuisance  parameter  is  the  actual  PDF  under  the  signal-absent  hypothesis. 

Application  of  the  technique  described  in  [1]  to  nonparametric  change  detection  first 
requires  assuming  that  the  change  (signal)  has  not  occurred  up  to  the  current  time  (referenced 
by  the  index  i).  A  sequential  test  between  the  null  hypothesis  (G(x)  =  F(x))  and  one  of 
the  above  alternatives  is  then  performed  where  data  distributed  under  F  are  obtained  from 
time  samples  i  -  K,  i  -  K  -  1,  . . . ,  i  -  K  -  L  +  I  where  K  represents  a  buffer  ensuring 
that  these  auxiliary  data  are  signal-free  and  L  is  the  maximum  number  of  available  auxiliary 
data  samples,  governed  by  the  stationarity  of  the  noise  sequence.  Data  distributed  under 
G  (which  is  either  F  under  the  signal-absent  hypothesis  or  one  of  the  above  alternatives  in 
the  signal-present  hypothesis)  are  then  obtained  from  time  samples  i-|-l,  i  +  2,  ...,  until 
a  decision  between  signal-absence  and  signal-presence  is  made,  say  at  time  index  i  +  M. 
If  a  signal- present  decision  is  made,  the  test  is  stopped  with  stopping  time  N  =  i  +  M. 
If  a  signal-absent  decision  is  made,  then  the  procedure  is  repeated  after  setting  i  =  i  + 
M.  Thus,  the  test  has  the  same  framework  as  the  Page  test  with  nuisance  parameter 
estimation,  and  results  in  the  same  ASN  of  equation  (39)  in  [1].  Evaluation  of  this  requires 
the  ASNs  and  operating  characteristic  of  the  individual  nonparametric  sequential  tests  used 
to  implement  this  nonparametric  change  point  detector,  which,  depending  on  the  test,  may 
not  be  straightforward  to  determine.  Important  in  this  evaluation  is  that  the  ASNs  and 
operating  characteristic  are  determined  accounting  for  the  randomness  of  the  auxiliary  data. 
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In  [1],  and  in  the  previous  sections  of  this  report,  this  is  done  by  conditioning  on  a  nuisance 
parameter  estimate  followed  by  an  expectation  removing  the  conditioning.  This  approach 
may  or  may  not  apply  to  any  given  nonparametric  sequential  test.  Govindarajulu  [9]  and 
Sen  in  [10]  describe  several  sequential  nonparametric  tests,  including  two-sample  tests  such 
as  that  of  Savage  and  Sethuraman. 

It  may  be  possible  to  modify  this  type  of  nonparametric  change  detection  problem  to  allow 
a  parameterized  noise  only  distribution  with  unknown  nuisance  parameters  and  sequential 
testing  for  some  sort  of  generic  change  indicating  signal-presence. 
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CONCLUSIONS 


The  Page  test  with  nuisance  parameter  estimation  [1]  was  applied  to  the  detection  of  the 
following  signal  types;  (i)  complex  Gaussian  distributed  data  with  a  shift  in  mean  signal 
and  unknown  variance,  (ii)  exponentially  distributed  data  with  an  unknown  noise  power, 
and  (iii)  multivariate  complex  Gaussian  distributed  data  with  a  deterministic  array  signal 
and  an  unknown  interference  covariance  matrix. 

For  each  signal  type,  a  detector  nonlinearity  was  developed  from  either  LO  or  GLR  prin¬ 
ciples  and  use  of  the  Dyson  bias  [4],  and  the  corresponding  Siegmund-based  approximations 
to  the  ASN  derived.  The  Siegmund-  and  Wald-based  approximations  were  compared  to  sim¬ 
ulation  for  the  average  time  between  false  alarms  and  average  delay  before  detection  where 
it  was  observed  that  the  Siegmund-based  approximations  are  more  accurate.  Corruption  of 
the  nuisance  parameter  estimate  by  auxiliary  data  containing  both  signal  and  noise  compo¬ 
nents  was  observed  for  weaker  signal  strengths.  This  may  be  avoided  by  proper  design  of  the 
buffer  zone  between  the  auxiliary  data  and  the  current  data  sample.  For  instance,  the  buffer 
size  may  be  chosen  so  that  minimal  corruption  occurs  for  the  smallest  signal  level  desired 
to  be  detected.  Operating  characteristic  curves  for  each  of  the  signal  types  were  generated 
describing  the  performance  of  the  detectors  as  a  function  of  threshold  and  SNR. 

The  probability  of  detecting  a  finite  duration  signal  was  approximated  by  the  Brownian 
motion  and  moment  matching  approximations  of  [2]  and  estimated  by  the  Poisson  mixture 
method  of  [3]  where  it  was  observed  that  the  Poisson  mixture  method  provided  the  minimum 
total  estimation  error.  This  is  an  expected  result  as  the  approximation  methods  of  [2]  are 
based  on  the  assumption  of  a  Page  test  update  sequence  that  is  independent  in  time  and,  in 
the  Brownian  motion  case,  has  Gaussian  update  statistics.  Use  of  the  nuisance  parameter 
estimate  in  the  Page  test  update  introduces  correlation  into  the  update  sequence.  The 
Poisson  mixture  method  simply  models  the  PDF  of  the  stopping  time  by  a  mixture  of 
Poisson  random  variables.  It  should  be  noted  that,  as  indicated  by  the  results  of  [11],  this 
approximation  may  be  improved  by  including  a  single  geometric  component  in  the  mixture 
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to  correctly  model  the  tail  probability. 

As  with  most  GLR  techniques,  the  only  optimality  associated  with  the  nonlinearity  choice 
for  the  above  signal  types  is  asymptotic  as  the  amount  of  auxiliary  data  becomes  infinite. 
Thus,  it  is  feasible  that  the  performance  of  these  signal  detectors  may  be  improved  by  a 
modification  of  the  nonlinearity.  For  instance,  the  asymptotically  optimal  bias  of  [4]  may 
be  used  rather  than  the  Dyson  bias.  However,  as  previously  mentioned,  this  results  in  a 
more  difficult  theoretical  analysis  of  the  detector  performance.  Other  modifications  may 
include  the  use  of  a  design  SNR  in  the  detector  implementation  rather  than  a  design  signal 
strength  or  inclusion  of  the  statistics  of  the  nuisance  parameter  estimate  into  the  nonlinearity 
formulation. 

Signal-onset  detection  where  signal-presence  is  distinguished  from  signal-absence  solely  by 
a  functional  relationship  between  the  CDFs  of  the  data  was  discussed  within  the  framework 
of  the  Page  test  with  nuisance  parameter  estimation. 
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APPENDIX  A:  MGF  UNITY  ROOT  FOR  A  NONCENTRAL 
CHI-SQUARED  RANDOM  VARIABLE 


Suppose  that  a  Page  test  update  has  the  form 

g{Y)  =  aY-h,  (A-1) 

where  a  >  0  and  the  data  Y  has  a  noncentral  chi-squared  distribution  with  n  degrees  of 
freedom  and  a  noncentrality  parameter  6. 

The  nonzero  unity  root  of  the  moment  generating  function  (MGF)  of  g  (Y)  is  the  value 
of  t  satisfying 


E  [( 


My) 


M{t)  =  l.  (A-2) 

The  MGF  of  g  may  be  obtained  from  the  MGF  of  a  noncentral  chi-squared  random  variable 

by 


M{t)  =  E 


j(aY-b) 


=  e"'’*E 

exp{-&^  +  lSt} 

(1  -  2at)^ 

iov  t  <  :A.  Define  (j)  (t)  as  the  natural  logarithm  of  M  (t), 


The  derivative  of  (t)  may  be  shown  to  be 


(A-3) 


(A-4) 


(A-5) 


(1  -  2atf 

It  may  also  be  shown  that  (j)"  (t)  >  0  for  alH  <  so  is  convex.  Clearly,  t  =  0  is  a 
root  of  4>  {t).  The  other  root  (t*),  which  is  the  nonzero  MGF  unity  root,  is  negative  (t*  <  0) 
when  the  mean  of  g, 


E[9(y)]  =  aE{y}-6 

=  a{n  +  S)  —  b  , 


(A-6) 

A-1 


TR  10,991 


is  positive  and  positive  {t*  >  0)  when  the  mean  is  negative.  Thus,  (p  (t)  may  look  like  the 
curves  found  in  figure  A-1. 

The  suggested  approach  to  finding  the  root  t*  is  to  find  the  minimum  point  of  0  (t),  move 
slightly  to  the  left  if  t*  <  0  or  slightly  toward  ^  if  t*  >  0  and  apply  the  Newton-Raphson 
iteration  for  root  finding  which  has  the  update  equation 


''fc+i 


0'W) 


(A-7) 


Termination  of  the  update  should  occur  when  successive  updates  of  t*  and  (pit*)  fail  to 
change  significantly. 

The  minimum  point  of  (p  (t)  occurs  when  cp'  (t)  =  0.  This  results  in  requiring 


t^  + 


b  —  a{n  +  6) 
^  4a% 


=  0 


(A-8) 


which  occurs  when 


t  =  t  1  lTV'l  +  (n-2  +  ^)^ 
^  2a  2b 


(A-9) 


If  n  >  2  or  n  =  1  and  5  >  1,  the  argument  of  the  square  root  in  equation  (A-9)  >  1.  Thus, 
since  valid  solutions  require  t  <  the  minimum  point  in  this  case  is 


t 


1  +  +  [n  -  2  +  S)  ^ 

— 


(A-10) 


For  the  case  of  n  =  1  and  0  <  5  <  1,  both  roots  should  be  computed  and  the  real  one  with 
the  correct  sign  chosen. 
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